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Abstract 

The inclusive t lepton hadronic decay is studied within Dispersive approach to QCD. The significance of effects due 
to hadronization is convincingly demonstrated. The approach on hand proves to be capable of describing experimental 
data on t lepton hadronic decay in vector and axial-vector channels. The vicinity of values of QCD scale parameter 
obtained in both channels bears witness to the self-consistency of developed approach. 
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The inclusive t lepton hadronic decay provides a 
clean environment for the study of the nonperturbative 
aspects of Quantum Chromodynamics (QCD) at low 
energies. In particular, this process is commonly em- 
ployed in tests of QCD and entire Standard Model, that, 
in turn, furnishes stringent constraints on possible New 
Physics beyond the latter. 

The relevant experimentally measurable quantity is 
the ratio of the total width of r lepton decay into hadrons 
to the width of its leptonic decay, which can be decom- 
posed into several parts: 

T(t~ — > hadrons - v r ) 

T = — 7Z— — 

= R% + R£ + R%+R£ + R„. (1) 

In the second line of this equation the first four terms 
account for the hadronic decay modes involving light 
quarks (u, d) only and associated with vector (V) and 
axial-vector (A) quark currents, respectively, whereas 
the last term accounts for the r lepton decay modes 
which involve strange quark. The superscript J indi- 
cates the angular momentum in the hadronic rest frame. 

All the quantities appearing in the second line of 
Eq. ([TJ can be evaluated by making use of the spectral 
functions, which are determined from the experiment. 
For the zero angular momentum (J = 0) the vector spec- 
tral function vanishes, whereas the axial-vector one is 
usually approximated by Dirac ^-function, since the 
main contribution comes from the pion pole here. The 
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experimental predictions [ 1 ] for the nonstrange spectral 
functions corresponding to J — 1 are presented in Fig.Q] 
In what follows we shall restrict ourselves to the consid- 
eration of terms R'° l w and R'° A of 7? T -ratio ([TJ. 

The theoretical prediction for the aforementioned 
quantities reads 

K% = y I Kdl 2 s EW (a^ a d + c) . (2) 

where N c = 3 is the number of colors, | V u &\ = 0.97425 + 
0.00022 is Cabibbo-Kobayashi-Maskawa matrix ele- 
ment |H,S E w = 1-0194 ±0.0050 and 6' m = 0.0010 stand 
for the electroweak corrections (see Refs. JHQl), and 

A^=2 f ' f(^-]R^( s )^- (3) 
M IA \Mf) M] 

denotes the QCD contribution. In this equation 

R v '\s) = - Im lim n v/A (i + is), (4) 

with YF IA (q 2 ) being the hadronic vacuum polarization 
function, and fix) — (1 - x) 2 (1 + 2x). 

In Eq. 0} Mi denotes the mass of the lepton on 
hand, whereas m v/A stands for the total mass of the 
lightest allowed hadronic decay mode of this lepton 
in the corresponding channel. The nonvanishing value 
of m v/A explicitly embodies the physical fact that 
t lepton is the only lepton which is heavy enough 
(M T =i 1.777 GeV 0) to decay into hadrons. Indeed, 
in the massless limit (m v/A = 0) the theoretical pre- 
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Figure 1: The inclusive vector and axial-vector spectral functions Q[. Vertical solid lines mark the boundaries of respective kinematic intervals, whereas horizontal 
dashed lines denote the naive massless parton model prediction. 



diction for A QCD © is nonvanishin^] for either lepton 
(/ = e, fi, t), whereas in the realistic case (m v/A + 0) 
Eq. Q acquires non-zero value for the case of the 
t lepton only. 

In general, it is convenient to perform the theoretical 
analysis of inclusive r lepton hadronic decay in terms 
of the Adler function [5] (the indices "V" and "A" will 
only be shown when relevant hereinafter) 



£>(fi 2 ) 



dTl(-Q 2 ) 
d\nQ 2 ' 



Q 1 



-9 



(5) 



Within perturbation theory the ultraviolet behavior of 
this function can be approximated by power series in 
the strong running coupling a s (Q 2 ): D(Q 2 ) D^ n (Q 2 ) 
for Q 2 — > co, where 



^ P 'en(G 2 ) 



1 + 



(6) 



1) «pen(fi 2 ) 



At the one-loop level (i.e., for I ■■ 
4n/(J3 lnz), z = Q 2 /A 2 ,/3 = 11- 2n f /3, A denotes the 
QCD scale parameter, «f is the number of active flavors, 
and d\ = l/n, see papers OaLZI] and references therein 
for the details. In what follows the one-loop level with 
rif = 3 active flavors will be assumed. 

For the beginning, let us study the massless limit, that 
implies that the masses of all final state particles are 
neglected (m = 0). By making use of definitions (|4|i 



1 Specifically, the leading-order term of Eq. (8) A"j~, = 1 (which 
corresponds to the naive massless parton model prediction for the 
Adler function {6) D^JQ 2 ) = 1) does not depend on M;, and, there- 
fore, is unique for either lepton. 



and ([5]), integrating by parts, and additionally employ- 
ing Cauchy theorem, the quantity A QCD (O can be repre- 
sented as (see Refs. [3, 8,|9[]) 

A QCD = ^Jd(m 2 e m ){\ +2e ie -2e m -e m )dG. (7) 

It is necessary to outline here that Eq. (0 can be de- 
rived from Eq. (O only for the massless limit of "gen- 
uine physical" Adler function D phys ((2 2 ), which pos- 
sesses the correct analytic properties in the kinematic 
variable Q 2 . However, in Eq. (0 one usually directly 
employs the perturbative approximation D p&n (Q 2 ) ©, 
which has unphysical singularities in Q 2 . At the 
one-loop level this prescription eventually leads to 



pert 



= A (0 U- 

peit A) 



4 f« AA\{6) + 0A 2 (0) 
n(A 2 + 2 ) 



d6, 



(8) 



where Ai(f9) = 1 + 2cos(f?) - 2cos(36») - cos(4f?), 
A 2 (0) = 2sin(<9) - 2sin(36») - sin(40), A = ln(M?/A 2 ), 

andAg, = l. 

It is worth noting also that perturbative approach pro- 
vides identical predictions for functions (0 in vector 
and axial-vector channels (i.e., Ap ert = Ap ert ). However, 
their experimental values Jll [loll are different, namely 



A^ xp = 1.224+0.050, 



A* xp = 0.748 +0.034. 



(9) 



The juxtaposition of these quantities with perturbative 
result ([8]) is presented in Fig. [2] As one can infer from 
this figure, for vector channel there are two solutions for 
the QCD scale parameter, namely, A = (434+|^)MeV 

and A = ( 1 652+23) Me V. Commonly, the first of these 
solutions is retained, whereas the second one is merely 
disregarded. As for the axial-vector channel, the 
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Figure 2: Comparison of the one-loop perturbative expression A peil {8} (solid curves) with relevant experimental data (9} (horizontal shaded bands). The solution for 
QCD scale parameter A (if exists) is shown by vertical dashed band. 



perturbative approach fails to describe the experimental 
data on inclusive t lepton hadronic decay, since for any 
value of A the function A pert (H) exceeds A* xp (|9}. 

It is crucial to emphasize that the presented above 
massless limit completely leaves out the effects due to 
hadronization, which play an important role in the stud- 
ies of the strong interaction processes at low energies. 
Specifically, the mathematical realization of the physi- 
cal fact, that in a strong interaction process no hadrons 
can be produced at energies below the total mass m of 
the lightest allowed hadronic final state, consists in the 
fact that the beginning of cut of corresponding hadronic 
vacuum polarization function Tl(q 2 ) in complex 
g 2 -plane is located at the threshold of hadronic produc- 
tion q 2 = m 2 , but not at q 2 = 0. Such restrictions are in- 
herently embodied within relevant dispersion relation. 
In turn, the latter imposes stringent physical nonper- 
turbative constraints on the quantities on hand, which 
should certainly be accounted for when one is trying to 
go beyond the limitations of perturbation theory. 

Thus, the nonperturbative constraints, which disper- 
sion relation [5] 



D(Q 2 ) 



Q- 



R(s) 



■ds 



(10) 



o + e 2 ) 2 

imposes on the Adler function (0, have been merged 
with corresponding perturbative result (0 in the frame- 
work of Dispersive approach to QCD, that has even- 
tually led to the following inte gral rep resentations for 
functions (|4]i and (0 (see Refs. 11 ill 1211 for the details): 



R(s) = r i0) (s) + 6\ 1 



dcr 



(11) 



f 00 A 

D(Q 2 ) = dP\Q 2 ) + P{Q\m 2 ,(f)p{ ( f) — . (12) 
Jm 2 °~ 



In these equations p(cr) denotes the spectral density, 
P(Q 2 ,m 2 ,cr) = Q 2 (o--m 2 )/((Q 2 +m 2 )(o-+Q 2 )), andfl(x) 
is the unit step-function (0(x) = 1 if x > and 9{x) = 
otherwise). It is worth noting also that in the massless 
limit (m = 0) expressions (fTTT i and ( fT2l become identi- 
cal to those of the so-called Analytic Perturbation The- 
ory |[l3l[l3l (see also Refs. lfl5l - [l7ll ). But, as it was men- 
tioned above, it is essential to keep the hadronic mass m 
nonvanishing within the approach on hand. 

Let us proceed now to the description of inclusive 
t lepton hadronic decay within Dispersive approach ifTll 
1211 (see also papers ESS and references therein). 
In this analysis the effects due to hadronization will be 
retained (in other words, the expressions ( TTTb and dT2l 
will be used instead of their perturbative approxima- 
tions and the hadronic mass m will be kept nonvan- 
ishing). The so-called "smooth kinematic threshold" 
for the leading-order term of function R{s) (see, e.g., 
Refs. II2UI22I1 ) will also be employed: 



(*) = (!- ml /*) 3/2 



<(e 2 )=i + 7 



l + ^ln 

2 



1+2^(1 -«(£)) 



(13) 
, (14) 



where m(^) = -\/l + and ^ = Q 2 /m 2 lk . Besides, the 
following model for the one-loop spectral density will 
be adopted: 



p(o-) = 



A 1 



A) ln 2 (o-/A 2 ) + n 2 cr 



(15) 



see papers 023142511 and references therein. The first term 
in the right-hand side of Eq. ( TT3T > is the one-loop pertur- 
bative contribution, whereas the second term represents 
intrinsically nonperturbative part of the spectral density. 



3 




0.00 0.25 0.50 0.75 0.00 0.25 0.50 0.75 



Figure 3: Comparison of expression Aq^ d 061 (solid curves) with relevant experimental data {9} (horizontal shaded bands). The solutions for QCD scale parameter A 
are shown by vertical dashed bands. 



Eventually, all this has led to the following expression 
for the quantity Aq C a d (01 within Dispersive approach (see 



Refs. [23, 26] for the details): 



A qcd = V 1 = £v/a (l + - ~€y IA + J^V/a) 



-3^ V /a| 1 + c <T V /a 22 ^ V/A 



x In 



M\ + vn^) - 1 

bV/A \ ' 

/ cr \ , , dcr 



(16) 



where H(x) = g(x)6(l-x)+g(l)e(x-l)-g({ v/A ), g(x) = 
x(2 - 2x 2 + x\ ml =* 0.075 GeV 2 , m 2 A 0.300 GeV 2 , 



and f Vj 



JM Z T - 



The comparison of obtained result (fTBT l with experi- 
mental data (|9]l yields nearly identical solutions for the 
QCD scale parameter A in both channels, see Fig. [3] 
Namely, A = (408 + 30) MeV for vector channel and 
A = (437 + 34) MeV for axial-vector one. Additionally, 
both these solutions agree with aforementioned pertur- 
bative solution for vector channel. 
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